Multigrid optimization for space-time discontinuous Galerkin discretizations of advection dominated flows

The goal of this research is to optimize multigrid methods for higher order accurate space-time discontinuous Galerkin discretizations. The main analysis tool is discrete Fourier analysis of two- and three-level multigrid algorithms. This gives the spectral radius of the error transformation operator which predicts the asymptotic rate of convergence of the multigrid algorithm. In the optimization process we therefore choose to minimize the spectral radius of the error transformation operator. We specifically consider optimizing h-multigrid methods with explicit Runge-Kutta type smoothers for second and third order accurate space-time discontinuous Galerkin finite element discretizations of the 2D advection-diffusion equation. The optimized schemes are compared with current h-multigrid techniques employing Runge-Kutta type smoothers. Also, the efficiency of h-, p- and hp-multigrid methods for solving the Euler equations of gas dynamics with a higher order accurate space-time DG method is investigated

On the design of block preconditioners for maritime engineering

The iterative error can be an important part of the total numerical error of any Com- putational Fluid Dynamics simulation when the iterative convergence stagnates or when loose convergence criteria are used. In the quest for better iterative convergence of CFD simulations, we consider the design of iterative methods for the Reynolds-averaged Navier-Stokes equations, discretized by finite-volume methods with cell-centered, co-located variables. The central point is the approximation of the Schur complement (pressure matrix) in the block factorization of the discrete system of mass and momentum equations. We show particular approximations of these blocks that yield either segregated solvers or block preconditioners for fully coupled solvers. The performance of these solvers are then demonstrated by computing the flow over a flat plate and around a tanker on both structured and unstructured grids. We find that iterative convergence to machine precision is attainable despite the high Reynolds numbers and mesh aspect ratio’s. Improved approximations of the Schur complement do result in improved convergence rates, but do not seem to pay-off in terms of total cost compared to the basic SIMPLE-type approximation

Extension of the discontinuous Galerkin finite element method to viscous rotor flow simulations

Heavy vibratory loading of rotorcraft is relevant for many operational aspects of helicopters, such as the structural life span of (rotating) components, op- erational availability, the pilot’s comfort, and the ef- fectiveness of weapon targeting systems. A precise understanding of the source of these vibrational loads has important consequences in these application ar- eas. Moreover, in order to exploit the full poten- tial offered by new vibration reduction technologies, current analysis tools need to be improved with re- spect to the level of physical modeling of flow phe- nomena which contribute to the vibratory loads. In this paper, a computational fluid dynamics tool for rotorcraft simulations based on first-principles flow physics is extended to enable the simulation of vis- cous flows. Viscous effects play a significant role in the aerodynamics of helicopter rotors in high-speed flight. The new model is applied to three-dimensional vortex flow and laminar dynamic stall. The applica- tions clearly demonstrate the capability of the new model to perform on deforming and adaptive meshes. This capability is essential for rotor simulations to accomodate the blade motions and to enhance vor- tex resolution

Space-time discontinuous Galerkin method for the compressible Navier-Stokes equations on deforming meshes

An overview is given of a space-time discontinuous Galerkin finite element method for the compressible Navier-Stokes equations. This method is well suited for problems with moving (free) boundaries which require the use of deforming elements. In addition, due to the local discretization, the space-time discontinuous Galerkin method is well suited for mesh adaptation and parallel computing. The algorithm is demonstrated with computations of the unsteady \ud ow field about a delta wing and a NACA0012 airfoil in rapid pitch up motion

Extension of a discontinuous Galerkin finite element method to viscous rotor flow simulations

Heavy vibratory loading of rotorcraft is relevant for many operational aspects of helicopters, such as the structural life span of (rotating) components, operational availability, the pilot's comfort, and the effectiveness of weapon targeting systems. A precise understanding of the source of these vibrational loads has important consequences in these application areas. Moreover, in order to exploit the full potential offered by new vibration reduction technologies, current analysis tools need to be improved with respect to the level of physical modeling of flow phenomena which contribute to the vibratory loads. In this paper, a computational fluid dynamics tool for rotorcraft simulations based on first-principles flow physics is extended to enable the simulation of viscous flows. Viscous effects play a significant role in the aerodynamics of helicopter rotors in high-speed flight. The new model is applied to three-dimensional vortex flow and laminar dynamic stall. The applications clearly demonstrate the capability of the new model to perform on deforming and adaptive meshes. This capability is essential for rotor simulations to accomodate the blade motions and to enhance vortex resolution

Chaotic multigrid methods for the solution of elliptic equations

Supercomputer power has been doubling approximately every 14 months for several decades, increasing the capabilities of scientific modelling at a similar rate. However, to utilize these machines effectively for applications such as computational fluid dynamics, improvements to strong scalability are required. Here, the particular focus is on semi-implicit, viscous-flow CFD, where the largest bottleneck to strong scalability is the parallel solution of the linear pressure-correction equation — an elliptic Poisson equation. State-of-the-art linear solvers, such as Krylov subspace or multigrid methods, provide excellent numerical performance for elliptic equations, but do not scale efficiently due to frequent synchronization between processes. Complete desynchronization is possible for basic, Jacobi-like solvers using the theory of ‘chaotic relaxations’. These non-deterministic, chaotic solvers scale superbly, as demonstrated herein, but lack the numerical performance to converge elliptic equations — even with the relatively lax convergence requirements of the example CFD application. However, these chaotic principles can also be applied to multigrid solvers. In this paper, a ‘chaotic-cycle’ algebraic multigrid method is described and implemented as an open-source library. It is tested on a model Poisson equation, and also within the context of CFD. Two CFD test cases are used: the canonical lid-driven cavity flow and the flow simulation of a ship (KVLCC2). The chaotic-cycle multigrid shows good scalability and numerical performance compared to classical V-, W- and F-cycles. On 2048 cores the chaotic-cycle multigrid solver performs up to faster than Flexible-GMRES and faster than classical V-cycle multigrid. Further improvements to chaotic-cycle multigrid can be made, relating to coarse-grid communications and desynchronized residual computations. It is expected that the chaotic-cycle multigrid could be applied to other scientific fields, wherever a scalable elliptic-equation solver is required

Space-time discontinuous Galerkin method for compressible flow

The space-time discontinuous Galerkin method allows the simulation of compressible flow in complex aerodynamical applications requiring moving, deforming and locally refined meshes. This thesis contains the space-time discretization of the physical model, a fully explicit solver for the resulting system of algebraic equations and all the necessary details for their implementation. The feasibility of the method is illustrated with various simulations, including the flow around a 3D delta wing and the dynamic stall of an airfoil in rapid pitch-up maneuvre

Runge-Kutta multigrid analysis for space-time discontinuous Galerkin discretization of an advection-diffusion equation

In this article, we analyse the convergence of multigrid (MG) iteration for solving the algebraic equations arising from a space-time discontinuous Galerkin (DG) discretization of the advection-diffusion equation. To keep the MG method fully explicit, we consider Runge-Kutta smoothers that solve the algebraic equations by marching in pseudo-time to steady state. Depending on the